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ABSTRACT 

We investigate the metal enrichment of elliptical galaxies in the framework of hierarchical 
models of galaxy formation. The semi-analytical model we use in this Letter, which has been 
used to study the metal enrichment of the intracluster medium (ICM) by Nagashima et al., in- 
cludes the effects of flows of gas and metals both into and out of galaxies and the processes of 
metal enrichment due to both type la and type II supernovae. We adopt a solar neighbourhood 
IMF for star formation in discs, but consider models in which starbursts have either a solar 
neighbourhood IMF or a top-heavy IMF. We find that the a-element abundance in ellipticals 
is consistent with observed values only if the top-heavy IMF is used. This result is consistent 
with our previous study on the metal enrichment of the ICM. We also discuss the abundance 
ratio of a elements to iron as a function of velocity dispersion and metallicity. We find that 
models with a top-heavy IMF match the a/Fe ratios observed in typical L„ ellipticals, but 
none of the models reproduce the observed increase of a/Fe with velocity dispersion. 
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1 INTRODUCTION 

The chemical abundances of elliptical galaxies have long been 
known to provide important constraints on theoretical models of 
their formation. The metallicities of the stellar populations in el- 
liptical galaxies are estimated observationally from the strength 
of stellar absorption features in their integrated spectra. Such 
measurements indicated an increase of the overall metal abun- 
dance with galaxy lumino sity and velocity dispersion (e.g. iFabeJ 
1 19731 iBender et alj Il993l) . and also an increase in the average 
Mg/Fe abundanc e ratio (e.g. lO'Connelll976tlWorthev et alJl992t 
J0rgensen 1999), with the central regions of most bright ellipti- 
cals having above-solar total metallicities and Mg/Fe ratios (the so- 
called a-enhancement). The detailed interpretation of line strengths 
in terms of chemical abundances is, however, complicated, due to 
the non-solar abundance ratios, and to the dependence of absorption 
line strengths on both metallicity and age. Detailed stellar popula- 
tion synthesis models have been developed to relate measured ab- 
sorption line strengths to ages and metallicities; older models (e.g. 
IWorthevI 1994b as sumed solar abundance rat ios, but newer models 
(e.g. lTraser et alj|2000fj: iThomas et alj|2003l) allow for non-solar 



ratios. Nonetheless, there remain significant uncertainties in going 
from measured line strengths to abundances of different elements. 

Most theoretical modelling of chemical abundances in ellipti- 
cals has been in the framework of the monolithic collapse or sin- 
gle burst model, in which all of the stars form at high redshift 
in a single burst of short duration, usually terminated by ejection 
of the remaining gas by a wind, following which the stars evolve 
passively (e.g. Larson 1975; Arimoto & Yoshii 1987). Subsequent 
works (e.g. lMatteucccilll994l : IThomas. Gre ggio & B ended [l 9991) 
have shown that it is possible to explain the trends of both increas- 
ing total metallicity and increasing Mg/Fe ratio with galaxy mass 
within this framework, provided that the star formation timescale 
of the initial burst (a free parameter) is chosen to have a suitable 
dependence on galaxy mass. 

However, there is now overwhelming evidence that struc- 
ture in the universe formed by hierarchical clustering of a 
dominant Cold Dark Matter (CDM) component, and strong 
evidence that galaxy mergers played an important role in 
the formation of elliptical galaxies. Galaxy formation within 
this hierarchical merger framework has been studied exten- 
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sively using semi-analytical models (e.g. White & Frenk 
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2002: Hatton et all2003l 


). There have, however, been very few de- 
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tailed studies of the ch emical evo lutio n of el liptical galaxies within 
this framework. [Cole et alj l200fjl) andfNagashim a & Yoshiil <2004 
computed the relation between total metallicity and luminosity for 
ellipticals, while Kauffmann & Chariot 1 1998) calculated the re- 
lation between Mg line strength and galaxy velocity dispersion, 
but these calculations all used the instantaneous recycling approx- 
imation for chemical evolution, with fixed solar abundance ratios 
built in. By construction, such models are unable to explain vari- 
ations in abund ance ratios such as Mg/Fe. Thomas 1 1999J) and 
iThomas & Kauffmanri Jl999l) calculated the separate evolution of 
[Fe/H] and [Mg/Fe], including both the prompt metal ejection by 
Type II supernovae (SNe II) and the long-delayed ejection by 
Type la supernovae (SNe la), based on star formation histories 
taken from a semi-analytical model. They found in their model that 
the average [Mg/Fe] actually decreased with increasing luminosity, 
in apparent conflict with observational data. However, they calcu- 
lated chemical evolution using a closed-box model, which ignored 
the effects of galaxy mergers and the transport of metals into and 
out of galaxies by gas inflows and outflows, all of which form part 
of standard semi-analytical models. 

In this paper we present the first calculations of chemical 
abundances in elliptical galaxies in a semi-analytical model which 
includes enrichment by both SNe la and SNe II within a fully con- 
sistent framework of hierarchical galaxy formation, i.e. including 
gas inflows and outflows and galaxy mergers. This enables us to in- 
vestigate whether hierarchical merger models are able to reproduce 
the observed variations in the ratio, a/Fe, of a-elements to iron. 
W e base our calculati ons on the semi-analytical model described 
in lBaugh et alj l2005t) . which assumes that stars formed in bursts 
triggered by galaxy mergers have a top-hea vy IMF , while stars 
formed quiescently in discs have a lKennicuttl ^83) IMF, similar 
to the solar neighbourhood. The top-heavy IMF in bursts was intro- 
duced in order tha t the model fit the number counts of faint sub-mm 
galaxies. Nagash ima et alj 12005) showed that this same model ex- 
plained the observed abundances in the intracluster medium (ICM) 
in galaxy clusters, when enrichment by SNe la was included. In the 
present paper, we compare predictions from the same mod el with 
observed abundances of elliptical galaxies. Nagashima & Okamoto 
1 2004) have already shown that a similar hierarchical model is able 
to explain the abundance patterns in disc galaxies like the Milky 
Way. 

The general idea that some or all of the stars in ellipti- 
cal galaxies may have formed with an IMF which was top- 
heavy rel ative to tha t in th e solar neighbourhood is n ot a 
new one. iMatteucccil ll994t). iGibson & Matteuccil 1 19971) and 
IThomas. Gre ggio & Benderl 7l999l) all proposed a top-heavv IMF 
as one possible explanation for the abundance patterns observed 
in ellipticals. The model presented here is novel in that it is based 
within a fully consistent hierarchical framework, and in that the pa- 
rameters of the top-heavy IMF and starbursts have been chosen in 
advance, independently of any metallicity constraints. We note that 
there is some more direct evidence for a top-heavy IMF in bursts 
from o bservations of the starburst galaxy M82 (Tsur u et alJll997t 
ISmith & GallagheJ200lh . 

Thus, the purpose of this Letter is to investigate theoretical 
predictions for the metallicities of elliptical galaxies from the same 
hierarchical galaxy form ation model as we used to study the en- 
richment of the ICM in Nagashima et al. (2005), and to make a 
preliminary comparison with observational data. In the next section 
we briefly summarise our model. We present the results from our 
model in Section 3 and our conclusions in Section 4. In the follow- 



ing, the cosmological parameters are fixed to be f2o = 0.3, Q,a = 
0.7, h = Ho/100 km s^Mpc -1 = 0.7 and cr 8 = 0.93. 



2 MODEL 

We adopt the same model as used in iNagashima et alj 120051) . 
which is an exte nsion of the GALFORM semi-analyt ical galaxy 
formation model <Cole et aljEoobl iBaugh et alj f2005) to include 
both prompt metal enrichment by SNe II and delayed enrichment 
by SNe la. The model, which will be described in full in Lacey 
et al. (2005, in preparation), includes the following physical pro- 
cesses: merging of dark haloes, radiative gas cooling, star forma- 
tion, SN feedback, mergers of galaxies due to dynamical friction, 
starbursts triggered by such mergers, and stellar population synthe- 
sis. It has been shown in our previous papers that this model repro- 
duces both the observed sub-mm source counts, and the observed 
metal abundances and baryon fractions in the ICM, pr ovided that 
a top-h eavy IMF is adopted for starbursts. We assume a K ennicuttl 
1 1983) IMF for quiescent star formation in galaxy discs, which has 
slope x = 0.4 for m < IMq and x = 1.5 for m > IMq. The 
top-heavy IMF has slope x — between 0.15 and 120Afo (where 
x = 1.35 for the Sa lpeter IMF). 

As shown by iBenson et alj 120031) . the original GALFORM 
model <Cole et alj |2000), when used with currently estimated 
cosmlogical parameters, predicted an excessive numbe r of ve ry lu- 
minous galaxies at the present epoch. Naga shima et alj|2005h con- 
sidered two alter native solutions to th is problem: (a) a superwind 
model [similar to lBenson et aljl2003l) 1. in which SN feedback ex- 
pels gas from haloes, the expelled gas later being recaptured once 
the halo grows up to a circular velocity V c = V raca ,p = 600 
km s _1 ; (b) a conduction model, in which the conductive heat 
flux from hot gas in the outer parts of massive haloes prevents gas 
from cooling in the centres of these haloes. Both models included 
a top-heavy IMF in bursts, and were consistent with metallicities 
in the ICM of clusters. We present results for both models in this 
paper. For comparison purposes, we also show results for a super- 
wind/Kennicutt model, which is identical to our standard superwind 
model except that a Kennicutt IMF is used throughout. 

We calculate the evolution of the abundances of different el- 
ements in the stars and gas using exactly the same prescription 
for chemical evolution as in Na gashima et alj l2005t) . The rates 
and chemical yields of SNe la and SNe II and gas restitution 
rates are computed consistently w ith t he a dopted IMFs ba sed o n 
IPortinari. Chiosi & BressarJl 19981) and lGreggio & Renzinilil983l) . 
We compute a total metallicity [Z/H] and iron abundance [Fe/H] 
(both being expressed as logarithmic values relative to solar). For 
comparing with the observational data, we also compute an a- 
element abundance [q/H] which includes O, Mg and other elements 
which are assum ed in observati onal analys es to va r y in s tep with 
these - seelTrager et al. 1 2000a) and Thomas et al. 1 2003) for de- 
tails. The a-abundance is dominated in mass fraction by O. The q- 
elements dominate the total metallicity, so [Z/H] is close to [a/H]. 

For comparing with observational data, we also need to com- 
pute stellar velocity dispersions of elliptical galaxies in our model. 
We compute the radii of stellar spheroids formed in mergers using 
the method described in lCole et alj l2000l) . and we then compute 
the circular velocity V c s at the half-mass radius of the spheroid in- 
cluding both its self-gravity and the gravity of its dark halo. We 
then estimate the ID stellar velocity dispersion as a = V c s / V3. 
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Table 1. Median of [a/H], [Fe/H] and [a/Fe] for galaxies with 150 ^ 
er/km s — 1 250. Values in the brackets represent the width between 25% 
and 75% values of the distributions. For the observations, we show both the 
measured central values and the estimated global values. The model values 
are all global values. 



model 


[a/H] 


[Fe/H] 


[a/Fe] 


observations (central) 


0.29 (0.19) 


0.06 (0.21) 


0.23 (0.10) 


observations (global) 


0.04(0.19) 


-0.19(0.21) 


0.23(0.10) 


superwind 


0.04(0.17) 


-0.18 (0.30) 


0.21 (0.14) 


superwind/Kennicutt 


-0.16(0.22) 


-0.27 (0.32) 


0.07(0.10) 


conduction 


0.29 (0.16) 


0.09 (0.28) 


0.21 (0.16) 



3 RESULTS 

Our predictions for the mean global stellar metallicities of ellipti- 
cals for the three different models we consider are shown in FigsQ 
and [2] We select ellipticals in the model as galaxies with B-band 
bulge-to-total light ratios (B/T)b > 0.6. The metallicities plot- 
ted for model galaxies are averages over all of the stars in the 
galaxy, weighted by the contribution of stars of different ages and 
metallicities to the V-band luminosity. This weighting is intended 
to mimic the measurement of metallicities from spectral features 
around 5000A. 

We have compiled a set of observational data to compare wit h 
the models, by combining the samples of [ Traeer e t alj |2000a). 
IProctor & Sansoml l20o3) , 1proctor et alJ <2004 and lThomas et alJ 
1 2005), and selecting ellipticals with measurement errors on [a/Fe] 
of less than 0.1 dex. This leaves us with about 140 ellipticals in a 
range of environments from groups to clusters, most of them with 
velocity dispersions in the range 10 < a < 300 kms -1 . (We 
note that a recent paper. iNelan et alJ 120051) . has reported very sim- 
ilar results to the above papers for much larger samples of 4097 
red-sequence galaxies in 93 low-redshift galaxy clusters.) In the 
above papers, metallicities have been derived from measured ab- 
sorption line strengths (the Lick indices) using stellar population 
models, assuming that all of the stars in a galaxy have identical ages 
and metallicities. The metallicities and velocity dispersions have 
been estimated from spectra of the central regions of the galaxies. 
Since there are radial line-strength gradients in elliptical galaxies, 
the central values differ from the global ave rage values. Based on 
the gradients measured bv lWu et alJ (2005), we estimate that the 
global total metallicities are typically lower by around 0.25 dex 
than the central values within 10% of the effective radius. In con- 
trast, the a/Fe ratio does not appear to show significant gradients 
(Meh lertet ai]|2003l) . Therefore, in comparing observational data 
with models, we include downward shifts by 0.25 dex in the ob- 
served values of [a/H] and [Fe/H], to correct them to global values, 
but assume that the observed central values of [a/Fe] are equal to 
the global values. 

The galaxies in the observational samples are selected in a 
heterogeneous way. In order to allow a fairer comparison between 
observations and models in FigsQand|2| we have plotted roughly 
the same number of model galaxies (~140) as in the observational 
sample, randomly selected from a volume-limited model catalogue 
so as to match the distribution of velocity dispersion a found in the 
observational sample. We note that the model reproduces the ob- 
served Lb-o correlation, where La is the total B-band luminosity. 

FigQshows [a/H], [Fe/H] and [a/Fe] plotted against velocity 
dispersion for the observational sample and for the three models 



we consider. The dashed boxes show the parameter region in which 
most of the observed elliptical galaxies are found, but in the panels 
showing the models, the boxes are plotted taking into account the 
estimated correction from central to global metallicities. To allow a 
more quantitative comparison, we also give in TableQthe median 
values and scatter for typical L» ellipticals, with a ~ 200km s~ . 
We note the following points: (1) The observed ellipticals show a 
large scatter in both abundances ([a/H] and [Fe/H]) and in abun- 
dance ratios ([a/Fe]) at any value of a, and this scatter is gener- 
ally reproduced by the models. (2) All three models predict simi- 
lar trends of increasing [a/H] with a, which agree with the trend 
seen in the observations. (3) However, the overall normalization of 
[a/H] differs between the models. When we take metallicity gra- 
dients into account, the superwind model (with a top-heavy IMF 
in bursts) seems to be in best agreement with the observed [a/H] 
(see Table 0. The superwind/Kennicutt model, which assumes a 
solar neighbourhood IMF for all stars, but is otherwise identical 
to the superwind model, predicts the [a/H] values too low by 0.2 
dex. The conduction model, which also has a top-heavy IMF in 
bursts, but has less gas and metals ejected from galaxies than in the 
superwind model, predicts [a/H] values too high by 0.2 dex. (4) 
The models all predict a trend of [Fe/H] increasing with a, in dis- 
agreement with the observational data, which indicate a flat depen- 
dence. (5) The superwind model predicts [Fe/H] for L» galaxies in 
agreement with observations, when metallicity gradients are taken 
into account (see Table 0. The superwind/Kennicutt model pre- 
dicts only slightly lower [Fe/H], showing that the top-heavy IMF 
in bursts affects the a-abundance more than the Fe-abundance. The 
conduction model predicts [Fe/H] too large by nearly 0.3 dex in L* 
galaxies. (6) The trends for [a/Fe] vs. a reflect those for [a/H] and 
[Fe/H] separately. The observations indicate a trend of [a/Fe] in- 
creasing with <7, while the models all predict a trend in the opposite 
sense. Both the superwind and conduction models predict super- 
solar [a/Fe] values for L* galaxies (a ~ 200 km s -1 ) which agree 
well with observations. However, none of the models match the 
typical observed [a/Fe] values in high-luminosity ellipticals, with 
a « 300km s" 1 . 

The trend we find in our models of the typical [a/Fe] de- 
creasing with increasing velocity dispe rsion or lumin osity is similar 
to what was found by Thomas & Kauffmann 1 1999). However, al- 
though |Thomas_& Kauffmann based their calculations on star for- 
mation histories extracted from hierarchical models, they computed 
the chemical evolution treating each galaxy as a closed box, ignor- 
ing galaxy mergers and inflows and outflows of gas and metals. 
Our calculation does include all of these effects, so our result is a 
stronger one. We will compare our results with the closed box re- 
sult in more detail in a future paper. In all three models presented 
here, ellipticals with a = 200 - 300 kms" 1 had around 40% of 
their star formation occuring in bursts, while the most recent burst 
typically happened 3-7 Gyr ago. 

In Fig|2| we plot the a/Fe ratio directly against the Fe- 
abundance, for exactly the same samples of observed and model 
galaxies as shown in FigQ We see that the observational data do 
not show any strong correlation of [a/Fe] with [Fe/H], but instead 
show a scatter of 0.4 dex in [a/Fe] at a given value of [Fe/H]. 
The models all show a trend of declining [a/Fe] with increasing 
[Fe/H], with varying amounts of scatter. Once we include the effect 
of metallicity gradients on the observational data, the superwind 
model seems to be the most consistent with observations, while the 
superwind/Kennicutt and conduction models predicts [a/Fe] values 
at a given [Fe/H] which are respectively too low or too high. 

We note here that there is a systematic uncertainty in the [a/H] 
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Figure 1. Abundances of a-elements, iron, and their ratio for stars in elliptical galaxies, plotted against velocity dispersion a. Abundances are on a logarithmic 
scale, relative to solar. The left-hand column shows observational data (not corrected for metallicity gradients) as follows: iThomasetaTlfcOOl T05, circles), 
IProctor et al J l2004l P04, sa»ares\ |Proctor & San som 1 2002, PS02, triangles), and lTrager et al J feoOOal TOO, crosses). The arrows show the estimated average 
corrections from central to global values, which are -0.25 dex for [a/H] and [Fe/H], and for [a/Fe]. For reference, dashed boxes surrounding most of the 
observational data are repeated across each row of panels. In the right 3 columns, these boxes are shifted by the estimated central-to-global corrections. The 
other 3 columns show predictions for the 3 different models, as indicated. The model abundances are V-band luminosity-weighted averages for each galaxy. 
The dotted boxes indicate how the observational data shift if the correction to [cr/Fe] proposed bv lProctor et al J i200 4 ) is adopted. 



and [a/Fe] values obtained from observational data due to tem- 
plate bias. The stellar population models used to derive metallici- 
ties and ages from the observed integrated absorption line strengths 
(Lick indices) are based on libraries of observed stellar spectra, in 
which some stars have non-solar a/Fe ratios. If no allowance is 
made for non-solar a/Fe in the template stellar spectra, then this 
can bias the determination of [a/Fe] and [a/H] from the integrated 
spectra of galaxies. However, different authors have made differ- 
ent corrections for this b ias, and the issue is currently controver- 
sial. Thomas et al. ( 2003) include an internal correction within their 
models. Pr octor et alj 12004) have instead proposed that the follow- 
ing correction should be applied to [a/Fe] values derived from in- 
tegrated line strengths: 

[a/Fe] flnal = -0.5[Fe/H] + [a/Fe] raw . (1) 

where [a/Fe] raw is the value derived using a stellar population 
model which ignores the template bias. A corresponding correc- 
tion is applied to [ a/H], but no correction is applied to [Fe/H]. 
IProctor et alj J200 4) found that when this correction is applied to 
observational data, the [a/H]-cr relation becomes flatter but tighter, 
while the trend of [q/Fe] inc reasing with a almost disappears. 
However, iThomas et alj J20051) h ave disputed t he validity of this 
correction and the conclusions of IProctor et all Since this correc- 
tion is controversial, we have plotted all of the observational data 



in Figs[T] and|2| (including the data from lProctor & Sanso nl <2002h 
and Procto ret alj 120041) ) without including this correction. How- 
ever, to show what the effect of making this correction would be, 
we have plotted in both figures dotted boxes indicating where the 
observational points would shift to if the correction were applied. 
We see that the model predictions in fact agree better with the cor- 
rected than uncorrected data, especially for the superwind model. 



4 CONCLUSIONS 

We have investigated the metal enrichment of elliptical galaxies in 
the framework of hierarchical galaxy formation models, taking into 
account the effects of galaxy mergers and gas inflows and outflows, 
as well as the production of metals by both SNe la and SNe II. This 
is the first time this has been done using a s emi-analytical model. 
Our model is the same as that which Naeashima et al. 1 2005j) used 
to show that the metallicities of iron and a elements in the ICM are 
successfully reproduced only when a top-heavy IMF is adopted for 
starbursts. We find that all of the models we consider predict corre- 
lations of the a-element abundance with stellar velocity dispersion 
a with similar slopes to what is observed, but the model without a 
top-heavy IMF in bursts (the superwinaVKennicutt model) predicts 
a-element abundances which are somewhat too low. Of our two 
models which include a top-heavy IMF in bursts, the superwind 
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Figure 2. The a/Fe abundance ratio vs. Fe-abundance for stars in elliptical 
galaxies, (a) Observational data, for the same galaxy samples as in Fig|2l 
Panels (b), (c) and (d) show predictions for the three different models. For 
reference, dashed boxes surrounding most of the observational data are plot- 
ted in all panels. In the panels showing the models, these boxes are shifted 
by the estimated central-to-global correction for metallicity gradients. The 
dotted parallelogram indicates how the observational data shift if the cor- 
rection to [a/Fe] proposed bv lProctor et alj 120041) is adopted. 



model predicts a-element abundances which agree well with ob- 
servations of bright elllipticals, while the conduction model (which 
has less ejection of gas and metals) predicts abundances which are 
too high. 

Our models with a top-heavy IMF also predict [a/Fe] values 
in typical typical L„ ellipticals (with a ~ 200 km s _1 ) which are 
similar to observed values, while the model in which all stars form 
with an IMF like that in the solar neighbourhood predicts [a/Fe] 
values which are too low. However, none of the models reproduce 
the trend of [a/Fe] increasing with a which is implied by most 
observational data. 

We note however some caveats to our comparison of the mod- 
els with observational data: (i) the observed metallicities are for the 
central regions of galaxies, while the models predict global val- 
ues; (ii) there are uncertainties in the estimation of metallicities 
from stellar absorption line indices connected with the treatment 
of non-solar a/Fe abundance ratios; (iii) the models used to esti- 
mate metallicities and ages from observed stellar absorption line 
indices may give biased results if a galaxy contains stars with a 
mixture of ages and metallicities, rather a single stellar population 
with a unique age and metallicity. We plan to address these issues 
in a future paper, to allow a more detailed test of our models against 
observational data on ellipticals. 



ACKNOWLEDGMENTS 

We acknowledge support from the PPARC rolling grant for extra- 
galactic astronomy and cosmology at Durham. MN and TO are sup- 
ported by the Japan Society for the Promotion of Science for Young 
Scientists (No.207 and 1891). CMB is supported by a Royal Soci- 
ety University Research Fellowship. We also thank the referee, Dr. 
Claudia Maraston, for a constructive and helpful report. 



REFERENCES 

Arimoto N., Yoshii Y, 1987, A&A, 173, 23 

Baugh CM., Lacey C.G., Frenk C.S. Granato GL. Silva L., Bres- 

san A., Benson A.J., Cole S., 2005, MNRAS, 356, 1191 
Bender, R.M., Burstein, D., Faber, S.M., 1993, ApJ 41 1, 153 
Benson A. J., Bower R.G, Frenk C. S., Lacey C. G, Baugh C. 

M., Cole S., 2003, ApJ, 599, 38 
Cole S., Aragon-Salamanca A., Frenk C. S., Navarro J. F., Zepf S. 

E„ 1994, MNRAS, 271,781 
Cole S., Lacey C. G, Baugh C. M., Frenk C. S„ 2000, MNRAS, 

319, 168 
Faber, S.M., 1973, ApJ 179, 731 
Gibson, B.K., Matteucci, F., 1997, MNRAS, 291, L8 
Greggio L., Renzini A., 1983, A&A, 118, 217 
Hatton, S., Devriendt, J.E.G, Ninin, S., Bouchet, F.R., Guider- 

doni, B., Vibert, D., 2003, MNRAS, 343, 75 
J0rgensen, I., 1999, MNRAS, 306, 607 
Kauffmann G, Chariot S., 1998, MNRAS, 294, 705 
Kauffmann G, White S. D. M., Guiderdoni B„ 1993, MNRAS, 

264, 201 

Kennicutt R.C., 1983, ApJ, 272, 54 
Larson, R.B., 1975, MNRAS 173, 671 
Matteucci, F, 1994, A&A, 288, 57 

Mehlert D., Thomas D., Saglia R.P., Bender R., Wegner G, 2003, 
A&A, 407, 423 

Menci N., Cavaliere A., Fontana A., Giallongo E., Poli F., 2002, 
ApJ, 575, 18 

Nagashima M., Lacey C.G., Baugh CM., Frenk C.S., Cole S., 

2005, MNRAS, 358, 1247 
Nagashima M., Okamoto T, 2004, submitted I astro-ph/0404486 1 
Nagashima M., Totani T, Gouda N., Yoshii Y, 2001, ApJ, 557, 

505 

Nagashima M., Yoshii Y, 2004, ApJ, 610, 23 

Nelan J.E., et al. 2005, MNRAS in press I astro-ph/0505301 1 

O'Connell, R.W., 1976, ApJ, 206, 730 

Portinari L., Chiosi C, Bressan A., 1998, A&A, 334, 505 

Proctor R.N., Forbes D.A., Hau G.K.T., Beasley M.A., De Silva 

G.M., Contreras R„ Terlevich AL, 2004, MNRAS, 349, 1381 
Proctor R.N., Sansom A.E., 2002, MNRAS, 333, 517 
Smith L.J., Gallagher III J.S., 2001, MNRAS, 326, 1027 
Somerville R.S., Primack J. R., 1999, MNRAS, 310, 1087 
Thomas D., 1999, MNRAS, 306, 655 

Thomas D., Greggio L., Bender R., 1999, MNRAS, 302, 537 
Thomas D., Kauffmann G, 1999, in Spectrophotometric dating of 

stars and galaxies, ed. I. Hubeny, S. Heap & R. Cornett, Vol. 192 

(ASP Conf. Ser.),261 
Thomas D., Maraston C, Bender R., 2003, MNRAS, 339, 897 
Thomas D., Maraston C, Bender R., De Oliveira CM., 2005, 

ApJ, 621, 637 

Trager S.C, Faber S.M., Worthey G, Jesus Gonzalez J., 2000a, 
ApJ, 119, 1645 

Trager S.C, Faber S.M., Worthey G, Jesus Gonzalez J., 2000b, 
ApJ, 120, 165 

Tsuru T.G., Awaki H., Koyama K., Ptak A., 1997, PASJ, 49, 619 
White S.D.M., Frenk C.S., 1991, ApJ, 379, 52 
Worthey, G, Faber, S.M., Jesus Gonzalez, J., 1992, ApJ 398, 69 
Worthey, G, 1994, ApJSupp, 95,107 

Wu H, Shao Z., Mo H.J., Xia X., Deng Z., 2005, ApJ, 622, 244 

This paper has been typeset from a TpX/ LTpX file prepared by the 
author. 



